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Abstract 

A theory for the nonlinear energy response of a system subjected to a heat bath is developed 
when the temperature of the heat bath is modulated sinusoidally. The theory is applied to a 
model glass forming system, where the landscape is assumed to have 20 basins and transition rates 
between basins obey a power law distribution. It is shown that the statistics of eigenvalues of the 
transition rate matrix, the glass transition temperature T g , the Vogel-Fulcher temperature Tq and 
the crossover temperature T x can be determined from the 1st- and 2nd-order ac specific heats, 
which are defined as coefficients of the 1st- and 2nd-order energy responses. The imaginary part of 
the lst-order ac specific heat has a broad peak corresponding to the distribution of the eigenvalues. 
When the temperature is decreased below T g , the frequency of the peak decreases and the width 
increases. Furthermore, the statistics of eigenvalues can be obtained from the frequency dependence 
of the lst-order ac specific heat. The 2nd-order ac specific heat shows extrema as a function of the 
frequency. The extrema diverge at the Vogel-Fulcher temperature To. The temperature dependence 
of the extrema changes significantly near T g and some extrema vanish near T x . 
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I. INTRODUCTION 



Since the anomaly of the specific heat at the glass transition was discovered in 1923 1 , 
many studies have been conducted to understand the behavior. Since the anomaly depends 
on the measurement process, the transition is now believed not to be understood in the 
framework of the standard thermodynamics^. 

Recently the concept of the landscape has been paid much attention because of possibility 
to explain the transition in non-equilibrium systems- 1 ^. In particular, the free energy land- 
scape (FEL) picture proposed by Odagaki et al 5,6 is considered to provide the unified concept 
for understanding thermodynamic and dynamic singularities of the glass transition. In fact, 
the single particle dynamics 7 - and the specific heat^^^^i^ were shown to phenomeno- 
logically be well described by the framework based on the FEL. Namely, the dynamical 
transition is understood as the Gaussian-to-Non- Gaussian transition 7 , and the thermody- 
namic singularities near T g , including the cooling rate dependence of the specific heat, is 
characterized as the Quenched-to- Annealed transition in dynamics on the FELr^i&ii 

The specific heat of non-equilibrium systems such as glass forming materials is defined by 
the reponse of the energy to the unit rise of temperature^. Because of the slow relaxations, 
the system cannot reach equilibrium during the measurement and the energy response shows 
the time delay. Therefore one can expect that the ac specific heat^^i 7 - will contain informa- 
tions of the slow dynamics. The ac specific heat Ci(iu>) is defined as the Laplace- Fourier 
transform of the energy or enthalpy correlation function <p(t}^2- 



The ac specific heat is often fitted by the Laplace- Fourier transform of the stretched exponen- 
tial 4>(t) = exp ^ — (t/r)^J, since the correlation function is believed to be the superposition 
of the Debye relaxation function with different relaxation times. However the origin of the 
distribution of the relaxation times has not clearly been understood. 

So far, studies of the energy response of glass forming systems have been limited within 
the linear response region, and the nonlinear energy response has not been studied yet. In a 
previous paper—, we proposed the description of the 1st- and 2nd-order energy response with 
the free energy landscape picture and applied it to a two level system with a diverging barrier. 
It is shown that the 2nd-order energy response has a diverging term at the temperature where 
the relaxation time diverges. 
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In this paper, we investigate the linear and nonlinear energy responses of non-equilibrium 
systems described by the FEL picture to an oscillating temperature and present the char- 
acteristic behavior of the 1st- and 2nd-order ac specific heat. Using a model FEL which 
supports a glass transition, we show that the Vogel-Fulcher temperature and the cross-over 
temperature as well as the glass transition temperature can be determined from the charac- 
teristic behavior of the ac specific heats. We also show that the statistics of the transition 
rate matrix representing the stochastic dynamics among the basins of the FEL can be ob- 
tained from the frequency dependence of the lst-order ac specific heat. 

The organization of the paper is as follows: In §2, we explain the 1st- and 2nd-order energy 
responses when the temperature of the heat bath is oscillated sinusoidally. The definition of 
the 1st- and 2nd-order ac specific heat is also given. In §3, the stochastic dynamics on the free 
energy landscape is explained. In §4, we describe the 1st- and 2nd-order energy responses 
to the oscillating temperature and the 1st- and 2nd-order ac specific heats when the system 
is described by the FEL picture. As an example, a model system with FEL consisting of 
20 basins is analyzed. In §5, where transition rates between basins are assumed to obey a 
power law distribution. We present the characteristics of the 1st- and 2nd-order ac specific 
heats of the glass former and show that the glass transition temperature, the Vogel-Fulcher 
temperature and the crossover temperature can be determined by the ac specific heats. In 
§6, our conclusion is given. 

II. THE FIRST AND SECOND ORDER SPECIFIC HEATS 

The specific heat at the constant volume Cy is conventionally defined by a derivative of 
the energy with respect to the temperature, 



where E cq is the energy in the equilibrium state, T is the temperature and V is the volume. 
In this discussion of the specific heat, it is not considered how long the system takes to 
equilibrate itself when the temperature is changed. When the system contains degrees of 
freedom of slow dynamics, one must consider the effect of the delay in response. We consider 
the energy response of a system with slow relaxations which is subjected to a heat bath whose 
temperature is oscillated as T + AT(t), where T is the average temperature and AT(t) is 
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the oscillating part. We assume that the energy response AE(t) can be expanded as follows: 



AE(t) = d^dit-t^ATit^ 

J —oo 

+ f dh [ rft 2 C 2 (t-t 1) t-t 2 )AT(t 1 )AT(t 2 ) + 0(AT 3 ), (3) 

J — oo J —oo 

where C% and C 2 represent the retardation effect of the system. 
The Fourier transform AE(u) of Eq. ([3]) is given by 

AE{u) = C , 1 (zcj)AT(cj) 

/oo 
duJxC 2 {iuJuiuJ - iuJi)AT{uj x )AT{uj - u x ) + 0(AT 3 ), (4) 
-oo 

where the Laplace components C\{p) and C 2 (pi,p 2 ) are defined by 

/■oo 

C 1 (p) = / dtC^e-P 1 , (5) 
Jo 

poo poo 

C 2 { Pl ,p 2 )= dh dt 2 C 2 {t x ,t 2 )e-^e-^\ (6) 
Jo Jo 

and AT(u) is the Fourier transform of AT(t). 

We now discuss the energy response AE(t) when AT(t) is a sinusoidal function AT(t) = 
T a sin(ujt), where T a is the amplitude of the oscillating temperature. It is straightforward to 
show that 

AE(t) = T a {C[{iu) sin(cut) + C'{(iu) cos(cut)} 
T 2 ~ 

- -^{C' 2 (iu, iu) cos(2cjt) - C' 2 \iu, iu) sin(2o;t) - C' 2 (iu, -iw)} + 0(T 3 ). (7) 

Here C\(t) and C 2 (ti,t 2 ) are assumed to be real and the notations ' and " represent the 
real and imaginary parts, respectively. The coefficients of the lst-order temperature term, 
C[(iu)) and C" (iu) are known as the real and imaginary parts of the (lst-order) ac specific 
heat, which was introduced by Birge and Nagel^ and by Christensen^. 

The 2nd-order temperature term consists of the oscillating and nonoscillating terms. We 
define the 2nd-order ac specific heat C' 2 {iuj,iuj) and C^ioo^iuo) by the coefficients of the 
oscillating term. 

The Laplace-Fourier transforms C\(p) and C 2 (p\,p 2 ) are related to the temperature 
derivative of the energy. When the change of the temperature is slower than that of the 
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time scale of the retardation effects C\(t) and 6*2(^1,^2), it is straightforward to obtain the 
following expressions, 

/°° BP 
dtC.it) = (8) 

/•OO /-OO 1 Cj2 Tf 

hm lim C 2 (p 1 ,p 2 )= / (ft / dhCifaU) = (9) 

When the frequency is smaller than the inverse of the structural relaxation time, the 
energy responses without delay and AE(t) can be expressed as 

op 1 a2 Tp 

AE(t) = ^AT(t) + ^^|?AT 2 (t) + 0(AT 3 ). (10) 

In this limit, C[(iu) becomes equal to the specific heat in the equilibrium state, Cy = 
8E eq /BT and C"(iu) vanishes. It can also be confirmed that in this limit C2{iuj,iu) and 
C^iu;, —iu) corresponds to the temperature derivative of the specific heat. These behaviors 
are consistent with Eqs. (jHJ) and (jHJ). 



III. DYNAMICS ON THE FREE ENERGY LANDSCAPE 

The free energy surface in the configurational space is defined by the partial partition 
function which is given by the partial summation of the phase space spanned by the fast 
microsopic motion 5 -^, and the slow dynamics is represented by the stochastic motion on 
the free energy surface. Around the glass transition point, the stochastic motions can be 
classified into two types; the fluctuation in one basin due to the structural fluctuation around 
a certain structure and the transition between basins which corresponds to the structural 
relaxation. 

Here, we concentrate on the transition between basins and ignore the fluctuation within a 
single basin. We denote the probability that the system is in basin i at time t at temperature 
T by Pi(T, t) and a physical quantity A of basin i by Ai(T). Then the physical quantity we 
measure at time t is defined as an average over the basins, 

A(T,t) = Y^Ai{T)Pi(T,t) = A(T) ■ P(T,t), (11) 

i 

where A(T) and P(T, t) are the vectors consisting of components Ai{T) and Pi(T, t), respec- 
tively. 
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The probability vector P(T, t) is assumed to obey the master equation, 

d 



dt 



P(T,t) = W(T)P(T,t), 



(12) 



where W(T) is the transition rate matrix, i.e. Wij is the transition rate from basin j to i 
and Wu = — Wji is the transition rate jumping out from basin i. The transition rate 
is assumed to be related to the free energy barrier as follows, 



W ij = Cexp[-P{F A (T)-F j {T))], 



(13) 



where C is the attempt frequency of the jump motion, (3 = l/ksT is the inverse of the 
temperature T multiplied by Boltzmann constant ks, Fi(T) is the free energy of basin i, 
Fa(T) = max{Fj(T), Fj(T)} + A^T) and Ajj(T) is the energy barrier between basin i 
and j. In the high temperature region, the transition rate is large and the system moves 
among basins freely. On the other hand, the transition rate becomes small and the structural 
transition is hindered at low temperatures. 

The transition rate matrix W in the master equation (fT2l) must satisfy the condition that 
the long time limit of Pi(T,t) becomes the Boltzmann distribution 

exph/^T)] 



}im P i (T,t) = P° q (T)- 
t-oo £\exp[-/3Fj(:r)] 

Therefore, the transition rate matrix W obeys the following detailed balance 



WP eq = 0. 



When the temperature T does not depend on time, Eq. (fl2"l) can be readily solved 



(14) 



(15) 



P(T,t) = exp[W(T)t]P(T,t = 0) 



V 



exp(Ait) 



\ 



V- x P{T,t = 0), 



(16) 



exp(AAr_it) 
^ exp(AArt) J 

where N is the number of basins, Aj is the z-th eigenvalue of W and V is a matrix whose 
columns are eigenvectors of W. It is important to note that there is an eigenvalue Aj = 
which corresponds to the detailed balance of Eq. (1T5T) . 
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IV. THE ENERGY RESPONSE TO THE OSCILLATING TEMPERATURE 



A. The 1st- and 2nd-order energy response 

The time dependent temperature T(t) is assumed to be T(t) = T + AT(t), where T is the 
average temperature and AT(t) is the oscillating part of the temperature. The probability 
vector P(T, t) is expanded as 

P{f{t),t) = P (t) + AP 1 {t) + AP 2 (t) + 0(AT 3 ), (17) 

where AP n is the term of order AT n . Note the explicit temerature dependence of the 
quantities on the right hand side of Eq. (TTTj) . is omitted. It is clear that Po(t) = P eq (T). 
Then the 1st- and 2nd-order energy responses, AEi(t) and AE 2 (t) at time t are represented 
as 

- dE - 

AE l (t) = E ■ AP l (t) + — • P eq (T)AT(t) (18) 

f) E 1 cP' E 

AE 2 (t) = E ■ AP 2 (t) + —AP^ATit) + ^P eq (T)AT 2 (t) (19) 

where E is the vector consisting of components E^. The 1st- and 2nd-order probability 
responses, AP.it) and AP 2 (t), obey the following equations derived from Eq. ([121) . 

jAP x {t) = W(T)AP 1 (t) + ^lp eq (T)AT(t) (20) 

d - dW(T) - 1 d 2 W(T) - 

-AP 2 (t) = W(T)AP 2 (t) + -^1 AP.it) AT {t) + --^2 Peq (T)AT(t) 3 . (21) 

It is straightforward to solve Eqs. (TSUI) and f[2~TT) as 

dse w ^ w P eq (T)AT(s), (22) 
AP 2 (t) = f ds x f l ds 2 e w «-^e w ^ ^P eq (T)AT( Sl )AT(s 2 ) 



' — oo J — oo 



+ U dse w ^ d -^AT\s). (23) 



B. The lst-order energy response and the ac specific heat 

The (lst-order) ac specific heat is defined as the linear coefficient of the energy response 
to the temperature change^. Comparing Eqs. (j3J), (JUJ) and ff22l . the retardation effect of 




FIG. 1: The schematic picture of the (lst-order) ac specific heat in the free energy landscape 
picture. 

the lst-order energy response C\(t) is represented in the free energy landscape picture as 
follows, 



d(t) = 2C (T)S(t) + E ■ exp(Wt) 



d_W p 



(24) 



where C (T) = P eq ■ dE/dT is the quenched specific heat contributed only from the fast 
degree of freedom at one basin^&i&iLM The 2nd term represents the response from the 
transition between basins. 

From Eq.(jS]), the lst-order ac specific heat C\{iuj) is given by the Laplace- Fourier trans- 
form of Ci(t). The real part C[(iu) and the imaginary part C"(iu) are written as 

\ 



( m<* + a?) 



C[{iu) = C {T) + E-V 
( 



\ 



dT 



(25) 



C'liiu) = E ■ V 



X 2 N /(^ + X%)J 
\ 



y- 



dT 



(26) 



y X n uj/(uj 2 + X%) J 

The frequency dependence of the ac specific heat is shown in FigfJJ schematically. In 
the low frequency region, C[(iuj) approaches the specific heat in the equilibrium state as 
C eq (T) = dE/dT -P eq + E- dP eq /dT and C'((iuo) vanishes. This indicates the fact that the 
energy responds without time delay when the temperature is oscillated slowly. Since the 
eigenvalue Aj is distributed, C"{iuj) has a broader peak than that of the Debye relaxation type 



and C[(iuj) is represented as the sum of Lorentzian funtions. This behavior is consistent with 
the measurement of the ac specific heat in glass forming materials^. In the high frequency 
region, C[(iu) approaches Cq(T) and C±(iuj) vanishes. It exhibits that the energy response 
within the lst-order temperature perturbation is determined by the energy change in a basin 
alone when the temperature fluctuates rapidly 

The lst-order ac specific heat can be used as the method to obtain the statistics of 
eigenvalues Aj. From Eqs. ff25]) and ff26]) . the scaled lst-order ac specific heat is related to 
the distribution of the eigenvalues A, of the transition rate matrix as follows, 

C eq -C -^ 9l \j + ^ {27} 

c^c Q -\ 9i ^^- (28) 

Here, <?j is a factor related to the vectors E, P eq and V given by 

gi = W* OT (29) 

<~<eq ~ ^0 

In the small frequency region, the imaginary part behaves as lim^Q C*i" (iuj)/(C eq — Co) = 
J2i giX[ l uj.\\i the high frequency region, the real part behaves as ]ixn u ^ 00 (C[(iu)) — Co)/(C eq — 
C ) = YliiQi'tfl^ 2 an d the imaginary part behaves as lim^oo Ci" {iuj)/(C eq — C ) = 
J^idi^i^ 1 - We can estimate Y^idi^ii Y^iQitf, an d Y^iQi^ 1 from the frequency depen- 
dence of the scaled lst-order ac specific heat. These properties correspond to (A), (A 2 ) and 
(A -1 ) respectively as will be shown in Sec. |Vj 

C. The 2nd-order ac specific heat 

The retardation effect of the 2nd-order energy response C 2 {ti, t 2 ) is represented from Eqs. 
(ED, (USD and ([23]) as 



*"1/T OTT7" OTT7" 

C 2 {t u t 2 ) = 2-^8{t 1 )8{t 2 ) + E ■ exp(H/t 1 )^rexp(H/t 2 )^ T P e ^(t 2 - t x ) 

1 - d 2 W - dE dW - 

+ -E ■ eMWh)—^^ - t 2 ) + 2— eMWh)—P^{t 2 ). (30) 

The 1st term expresses the response due to the instant change of the specific heat of each 
basin and the remaining terms correspond to the effect of the transition motion between 
basins. Here, 5(x) is the Dirac's delta function and 6(x) is the Heaviside step function. 



From the definition ([6]), the Laplace transform C^Qoi, P2) is represented as 
1 c)C! -> „ f)W * r)P eiJ 

C 2 (P1,P2) = - P2 B ' ((P, + P,)l - W)" l ^-(P2l " ^)"'^ 

_ . +p2)i _ W) -^ _ pi | fei _ (31) 

The 2nd-order ac specific heats, which are introduced in eq. ([3]) by the coefficients of the 
oscillating term in the 2nd-order energy response, C^^oj, ^) = C' 2 (iu, ^0 + iC 2 (iuj, iuj) are 
expressed as 

C' 2 (iu>, iu) = - 2co 2 E ■ (W 2 + ^Y x -^r 



+ • (W 2 + Au 2 )- l {2^-W + W^-\(W 2 + lu 2 Y 



dT 



-^ir (w V' (32) 

C' 2 \ito, iu) = uE ■ W(W 2 + Auj2 )~ 1 -gf 2 - 

+ a ,g (W ^ + ^)-i^l. (33) 

In the low frequency limit, C' 2 (iuj,iuj) becomes (dC eq /dT)/2 and C 2 (iu,iu) vanishes. In 
the high frequency limit, C 2 (iu>,iuj) is equal to (dC /dT)/2 and C 2 (iu,iu) vanishes. These 
behaviors are qualitatively similar to those of the lst-order ac specific heat. It indicates that 
the structural change of the system can catch up with the temperature change in the small 
frequency region and can not in the high frequency region. 

Terms including the differential of the transition rate with respect to the temperature, 
dW/ dT, is expected to show the characteristic behavior in the low temperature region, where 
the transition matrix W depends strongly on the temperature change. This characteristic 
does not appear in the lst-order ac specific heat. In the high temperature region, the tem- 
perature dependence of the transition rate matrix W is weak. Therefore, the terms including 
dW/dT are negligible, and the 2nd-order ac specific heat becomes the superposition of the 
susceptibility of the Debye relaxation. We explain this behavior for a trap model in the next 
section. 

The 2nd-order energy response can be measured in the materials, where the derivatives 
of the quenched and equilibrium specific heats are much different, since the order of the 2nd- 
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order energy response is the difference between derivatives of these specific heats, d(C eq — 



V. APPLICATION TO A MODEL FREE ENERGY LANDSCAPE OF GLASS 
FORMING SYSTEMS 

A. Model 

To clarify the characteristics of the 1st- and 2nd-order ac specific heat of glass forming 
materials, we apply the present analysis of the ac specific heat to a model landscape in- 
corporated with the distribution of the transition rate exploited in the trapping diffusion 



We prepare a landscape consisting of 20 basins which are mutually connected. The 
microscopic dynamics of the system in each basin is assumed to be the Debye oscillator. 
The free energy, F±, of basin i is given by 



where the energy e, L of the bottom of basin i is uniformly distributed with the variance e 
and Qjj is the Debye temperature. Note that the form of the distribution does not play any 
important roles in the following discussion. 

The energy, E{, of basin % is calculated from the temperature derivative of the free energy 
Fi. The transition rate w is assumed to obey the power law distribution P(w) as in the 
trapping diffusion modet^, 



C )/dT. 



model^ 1 ^. 




(34) 




otherwise 




(35) 



where p is related to the configurational entropy s c (T) as 



P 



TS C (T) - T g S c (T g ) 
T 9 S c (T g ) 



(36) 



The transition rate from basin j to i is now a random variable given by 



Wij = w x—i exp[-P{F A (T) - Fj(T))]. 



(37) 
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FIG. 2: The scaled ac specific heat in the trap model at the temperature T/T g = 0.85 (solid line), 
1.0 (dot line), 1.25 (dash line): 

Here, x is the uniform random number in [0, 1]. 

For numerical calcuations in this study, the exponent p of the jump rate distribution is 
simplified as 

T — T 

P=±—±, (38) 

i.e. S C (T) is assumed as TS C (T) ocT — T with the Vogel-Fulcher temperature T ~ Tk and 
the parameters are chosen as T /T g = 0.75 and ®r>/T g = 12.5, respectively. In addition, the 
crossover temperature T x is T x /T g = 1.25 since the crossover temperature is identified with 
the exponent p = 1.0^. 



B. The lst-order ac specific heat 

The real and imaginary parts of the scaled lst-order ac specific heat (C\{iu) — Co) / {C eq — 
Co) are shown in FigO As the temperature is reduced, a peak of the imaginary part shifts 
to the low frequency region and the width of the peak increases. 

Figure [3] shows the temperature dependence of the frequency oJ pea k at the peak of the 
imaginary part. In the low temperature region below T g , oj pea k obeys the power law func- 
tion as <jJ ve ak ~ (T — T ) a914 . In the high temperature region above T g , the temperature 

12 



o 




(T-T )/T g 

FIG. 3: The temperature dependence of the frequency at the peak of the imaginary part of the lst- 
order ac specific heat: The peak frequency (black circle) becomes small as the temperature reduces 
toward the Vogel-Fulcher temperature To. The line is fitted with the data in the low temperature 
region. 

dependence changes and u) pea k approaches toward 20w . It is straightforward to show that, 
when N basins are mutually connected, the uj pea k becomes Nwq in the high temperature 
region. In fact from Eq. fl37j) . all transition rates are wq when the temperature is much 
higher than the glass transition temperature T g . The relaxation time due to eigenvalues 
of transition matrix W is not distributed and the lst-order ac specific heat becomes the 
same as the susceptibility of the Debye relaxation. The master equation (|T2l) can be solved 
without difficulties in this limit as P^T, t) = e- Nwot (Pi(t = 0) - l/N) + l/N. Thus the peak 
frequency ui pea k becomes Nw in the high temperature region. 

We define the stretching factor a by the ratio a = u + /ui- (u + > uJ). Here, u± is the 
frequency satisfying the relation C"(iu±) = Ci(iu pea k)/2. Figure H] shows the temperature 
dependence of the stretching factor a as a function of temperature. In the high temperature 
region above T g , a is the same as that of the Debye susceptibility. As the temperature is 
reduced toward T , a of the peak increases and diverges at T . It means that the relaxation 
time is distributed and the broad peak appears in C"{iu) below T g . It is interesting to note 
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FIG. 4: The temperature dependence of the stretching factor a of the imaginary part. The line 
represents the fitted line by Vogel-Fulcher law. Here, u Debye = 

(2 + \/3)/(2 - \/3) is the stretching 
factor for the Debye relaxation type. As the temperature is reduced toward the Vogel-Fulcher 
temperature To, a increases. 

that a can be fitted well by the Vogel-Fulcher law as a ~ exp(A/(T — T )) shown in FigJH 
The statistics of eigenvalues of the transition rate matrix W can be obtained 
from the frequency dependence of Ci(za;).(See §IVBI ) We show a result of com- 
parison of \im UJ -> 00 ujC" 1 (iu)/(C eq - C ), lim^^u 2 (C[(iu) - C )/(C eq - C ), and 
lim^o cu^C" i(iu>) / (C eg — C ) with statistics of eigenvalues (A), (A 2 ) and (A -1 ) in FigfSl 
Here, (/(A)) represents the statistics of eigenvalues as (/(A)) = Ylf=i x^o /(-^)/20, which is 
calculated numerically from the transition matrix W. The result indicates that the contri- 
bution of gt is negligible and the coefficients obtained from the frequency dependence of the 
lst-order ac specific heat can be used to estimate the statistics of eigenvalues of transition 
rate matrix. 

Figure M shows the temperature dependence of lim^oo ljC" i{iuj)/{C eq — C )/W , 
hm^ooK^M - C )/(C eq - C )/W^ 2 , and lim^ ^(C e , - C ) / (C[(icu) - C )/W 
and the deviation between these statistics becomes significant below T g . It reflects that the 
distribution of the eigenvalues becomes broader below T g . 
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(T-T )/T g (T-T )/T g (T-T )/T g 

(a) (b) (c) 

FIG. 5: The temperature dependence of (a)— Hindoo uC" / (C eq — C$)/Wq, 

(bJIim^^ooW^CKiw) - C )/(C eq - C )/W 2 , and (c)-Wblim w -»ow- 1 (C w i(ia;) - C )/(C eq - C„) 
determined from the frequency dependence of the lst-order ac specific heat. The solid lines 
represent (a) — (A) /Wo, (d)(A 2 )/Wq and (c)— Wo(A -1 ), which are calculated numerically from 
eigenvalues of the transition matrix W, respectively. 




FIG. 6: The temperature dependence of — lim^^oo u;C"i(iLj)/(C eg — Co)/Wo, lira UJ ^ 00 {oj 2 (C' 1 (ioj) — 
C )/(C eq - C )/W 2 } 1 / 2 , and - lim^ uj(C eq - C )/{C[{iu) - C )/W . The deviation of these 
properties becomes significant below T g . 
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FIG. 7: The real part of the scaled 2nd-order ac specific heat at temperatures T/T g =0.825, 0.875, 
1.0, 1.25 and 2.5: As the temperature is decreased below the temperature Tx/T g = 1.25, a local 
minimum decreases. 

C. The 2nd-order ac specific heat 

The real and imaginary parts of the scaled 2nd-order ac specific heat {C2(ioj,iui) — 
€2(100, ioo))/(C , 2(0, 0) — 62(200, ioo)) are shown in Fig. [7] and Fig. [HI respectively. The 
real part shows a decreasing local minimum when the temperature is reduced below the 
crossover temperature Tx/T g = 1.25. It is important to note that the similar behavior is 
shown for the imaginary part in FigJHJ where the imaginary part has a decreasing minimum 
and an increasing maximum below the crossover temperature. 

The deviation from the Debye relaxation appears clearly in the Cole-Cole plot of the 
2nd-order ac specific heat shown in FigJHl In the high temperature region, the plots show 
the semicircle, since the 2nd-order ac specific heat is equal to the susceptibility of the Debye 
relaxation. As the temperature is reduced, the plots deviate from the semicircle, that shows 
that the terms including dW j dT in Eqs. (1521 and (1551) give significant contributions. 

The temperature dependence of the extrema of the 2nd-order ac specific heat is shown 
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FIG. 8: The imaginary part of the scaled 2nd-order ac specific heat at temperatures T/T g =0.825, 
0.875, 1.0, 1.25 and 2.5: As the temperature is decreased below the temperature Tx/T g = 1.25, a 
local minimum decreases and a local maximum increases. 

in FigrjJJ Below T g , the extrema obey the power law as (T — T )~ u , as the temperature is 
decreased toward To. The exponents v are v = 0.97, 1.05, 1.01 for the local maximum of the 
real part, the local minimum and the local maximum of the imaginary part, respectively. 
Above T g , the temperature dependence of the minimum of the imaginary part deviates from 
the power law and approaches to a constant value. 

VI. CONCLUSION 

We presented the theoretical formalism of calculating the energy response to the oscil- 
lating temperature in the free energy landscape picture. The characteristic behaviors of the 
energy responses are summarized as follows. 

The large and small frequency limits of the lst-order ac specific heat are identical to the 
equilibrium and quenched specific heat in the free energy landscape picture, respectively. 
The distribution of the eigenvalues of the transition rate matrix among basins is shown to 
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FIG. 9: The Cole-Cole plot of the scaled 2nd-order ac specific heat at temperatures T/T g =0.875, 
1.0, 1.25 and 2.5: At a high temperature, the plots are semicircle. As the temperature is decreased, 
the plots deviate from the semicircle. 

give rise to the broad peak of the imaginary part of the lst-order ac specific heat. This is 
consistent with the experiment for the glass forming materials. 

The 2nd-order ac specific heat is defined for the first time as the coefficient of the oscil- 
lating part of the 2nd-order energy response. The large and small frequency limits of the 
2nd-order ac specific heat are identical to the temperature derivative of the annealed and 
quenched specific heat, respectively. 

We analyzed the 1st- and 2nd-order ac specific heats for a model glass former on the basis 
of the free energy landscape picture and showed characteristics of the ac specific heats. 

A frequency uj pea k at a peak of an imaginary part of the lst-order ac specific heat shifts 
to the low frequency region as the temperature is decreased below the glass temperature 
temperature T g and approaches to at the Vogel-Fulcher temperature T . In the high 
temperature region above T g , u pea k becomes virtually independent of the temperature and 
approaches a constant value. A width of the peak becomes larger than that of the Debye 
relaxation type below T g due to the distribution of the eigenvalue A and diverges at To. The 
statistics of A such as (A), (A 2 ), and (A -1 ) can be obtained from the frequency dependence 
of the lst-order ac specific heat, and these quantities shows deviation from the Debye-type 
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FIG. 10: The temperature dependence of the extrema of the scaled 2nd-order ac specific heat: 
The lines are fitted by power law functions in the low temperature region. Below T g , the extrema 
obeys the power law as (T — Tq)~ u '. Above T g , the temperature dependence of the extrema shows 
significant changes. 

relaxation below T g . 

The scaled 2nd-order ac specific heats have extrema, which obey the power law as (T — 
Tq) u below T g . Above T g , the temperature dependence of the extrema changes. The local 
minimum of the real part and the local maximum of the imaginary part vanish and the 
local minimum of the imaginary part approaches to the constant value near the crossover 
temperature T x . 

In conclusion, the characteristics of the transition rate among basins reflect the frequency 
dependence of the 1st- and 2nd-order ac specific heats. It means that the glass transition 
point T g , the Vogel-Fulcher temperature T and the crossover temperature T x can be de- 
termined from the 1st- and 2nd-order ac specific heats. It is interesting to note that we can 
determine T , at which the shear viscosity diverges, i.e. a dynamically denned character- 
istic temperature, from the thermodynamic measurement as the 1st- and 2nd-order energy 
response. 
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